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c/3 ' A set of reduced, 2-D, two-fluid, drift-MHD (magnetohydrodynami- 

C^ ' 

0^' cal) equations is derived. Using these equations, a complete and fully self- 

O ■ consistent solution is obtained for an isolated magnetic island propagating 

^~»' through a slab plasma with uniform but different ion and electron fluid veloc- 

^ • 
P i' 

ities. The ion and electron fluid flow profiles around the island are uniquely de- 

K^ termined, and are everywhere continuous. Moreover, the island phase- velocity 

(N 

fr^ ■ is uniquely specified by the condition that there be zero net electromagnetic 

^H \ 

^«0 ■ force acting on the island. Finally, the ion polarization current correction to 



o 



X 



the Rutherford island width evolution equation is evaluated, and found to be 



Q ■ stabilizing provided that the anomalous perpendicular ion viscosity signifi- 

^' cantly exceeds the anomalous perpendicular electron viscosity. 



I. INTRODUCTION 

Tearing modes are magnetohydrodynamical (MHD) instabilities which often limit fusion 
plasma performance in magnetic confinement devices relying on nested toroidal magnetic 
flux-surfaces.^ As the name suggests, "tearing" modes tear and reconnect magnetic field- 
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lines, in the process converting nested toroidal flux-surfaces into helical magnetic islands. 
Such islands degrade plasma confinement because heat and particles are able to travel ra- 
dially from one side of an island to another by flowing along magnetic field-lines, which is a 
relatively fast process, instead of having to diffuse across magnetic flux-surfaces, which is a 
relatively slow process.^ 

Magnetic island physics is very well understood within the context of single-fluid MHD 
theory. According to this theory, the island width is governed by the well-known nonlinear 
evolution equation due to Rutherford.^ Moreover, the island is required to propagate at the 
local flow velocity of the MHD fluid, since fluid flow across the island separatrix is effectively 
prohibited. 

Magnetic island physics is less completely understood within the context of two-fluid, 
drift-MHD theory, ^^^^ which is far more relevant to present-day magnetic confinement de- 
vices than single-fluid theory. In two-fluid theory, the island is generally embedded within 
ion and electron fluids which flow at different velocities. The island itself usually propagates 
at some intermediate velocity. For sufficiently wide islands, both fluids are required to flow 
at the island propagation velocity in the region lying within the island separatrix (since 
neither fluid can easily cross the separatrix). However, the region immediately outside the 
separatrix is characterized by strongly sheared ion and electron fluid flow proflles, as the 
velocities of both fluids adjust to their unperturbed values far away from the island. The 
polarization current generated by the strongly sheared ion flow around the island separatrix 
gives rise to an additional term in the Rutherford island width evolution equation, which 
is stabilizing or destabilizing, depending on the island propagation velocity relative to the 
unperturbed flow velocities of the ion and MHD fluids. The key problems in two-fluid island 
theory are the unambiguous determination of the island phase-velocity, and the calculation 
of the ion and electron fluid flow proflles around the island separatrix. As yet, no consensus 
has emerged within the magnetic fusion community regarding the solution of these problems. 

In this paper, we flrst develop a set of reduced, 2-D, two-fluid, drift-MHD equations. 
These equations contain both electron and ion diamagnetic effects (including the contri- 



bution of the ion gyroviscous tensor), as well as the Hall effect and parallel electron com- 
pressibility. However, they do not contain electron inertia or the compressible Alfven wave 
(which play negligible roles in conventional magnetic island physics). Our set of equations 
consist of four coupled partial differential equations, and is both analytically tractable and 
easy to solve numerically. We employ our equations to study the evolution of an isolated 
magnetic island in slab geometry. Using a particular ordering scheme, we are able to cal- 
culate the island phase-velocity, and to uniquely determine the ion and electron fluid flow 
profiles outside the island separatrix. 

II. DERIVATION OF REDUCED EQUATIONS 

A. Introduction 

In this section, we shall generalize the analysis of Refs. 19 and 20 to obtain a set of 
reduced, 2-D, two-fluid, drift-MHD equations which take ion diamagnetic flows into account. 

B. Basic equations 

Standard right-handed Cartesian coordinates (x, y, z) are adopted. Consider a quasi- 
neutral plasma with singly-charged ions. The ion/electron number density no is assumed to 
be uniform and constant. Suppose that Tj = rTe, where Tj^g is the ion/electron temperature, 
and r is uniform and constant. 

Broadly following Ref. 21, we adopt the following set of two-fiuid, drift-MHD equations: 
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^ + V-v]p = -TPV-V + kV^P. (3) 



Here, E is the electric field, B the magnetic field, J the electric current density, V the plasma 
guiding-center velocity, P the total plasma pressure, e the magnitude of the electron charge, 
rrii the ion mass, rj the (uniform) plasma resistivity, fig the (uniform) electron viscosity, 
fii the (uniform) ion viscosity, k, the (uniform) plasma thermal conductivity, and F = 5/3 
the plasma ratio of specific heats. Furthermore, b = B/B, and V^ = b x VP/erioB. 
The above equations take into account the anisotropic ion gyroviscous tensor, but neglect 
electron inertia. Our system of equations is completed by Maxwell's equations: V-JB = 0, 
V X E = —dB/dt, and V x i? = /iq J. Note that the transport coefficients, /ij, fie, and 
K,, appearing in the above equations, are phenomenological in nature, and are supposed 
to represent the anomalous diffusive transport of energy and momentum across magnetic 
fiux-surfaces due to small-scale plasma turbulence 

C. Normalized equations 

Let V = aV, t = t/(a/K), B = B/B^, E = E/{BaVa), J = J/{Ba/fioa), V = 
V/Va, P = P/{B^/fio), V = v/ifJ'oVaa), fii,e = fJ,i,e/ {noniiVa a) , k = K/iVaa), where 
Va = Ba/ ^HonQirii. Here, a is a convenient scale length, and B^ a convenient scale magnetic 
field-strength. 

Neglecting hats, our normalized two-fiuid equations take the form: 

E + VxB + di(vP-Y^ib-\/P)b-JxB-iXe^We) =^(''- 2TT7^*)' ^^'' 

f d „ diT \ diT „ , , , , „ ^ 

TST + ^-V + -r^— K-V y - -^— V,-V{[b-V] b)=JxB-VP 
\dt l + T J l + T 

+/iiVV, + /ieWe, (5) 

I^ + V-v]p = -TPV-V + kV^P. (6) 

Here, V^ = b x VP/B, and di = (mj/no e^/io)^/^/a is the normalized collisionless ion skin- 
depth. Maxwell's equations are written: V-B = 0, V x E = —dB/dt, and V x B = J. 



D. 2-D assumption 

Let us make the simplifying assumption that there is no variation of quantities in the z- 
direction: i.e., djdz = 0. It immediately follows that B = Vip x z + Bz z, and Ez = —dip/dt. 

E. Reduction process 

Let us adopt the following ordering, which is designed to decouple the compressional 
Alfven wave from all the other waves in the system: 

P = Po + 5oPi+P2, (7) 

Bz = Bo + bz. (8) 

Here, Pq and Bq are uniform and constant, and 

Po > 5o > 1. (9) 

Furthermore, pi, p2, bz, ip, V, V, and d/dt are all assumed to be 0(1), and V-V is assumed 
to be much less than 0(1). 

Now, to lowest order, the 2;-component of Ohm's law, Eq. (4), gives 

[^ + V- V U = -d, [bz + rpi/(l + r),^ +r]V^ij- d, /ie V\Vz + d, V^)- (10) 

Here, [^, P] = VA x VB-z. Likewise, the z-component of the curl of Eq. (4) reduces to 

(^ + V-VJ6, = [Vz + d,V^ij,ij] - BoV -V + r]V^(^bz + ^Y^Pi 



+di fie V^ 



U-diV'(bz + -^—Pi 
1 + r 



(11) 



Here, U = -V x V-z. 

To lowest order, the equation of motion, Eq. (5), implies that 

Pi ^ -bz. (12) 

Furthermore, the ^-component of this equation yields 



dt 
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whereas the z-component of its curl reduces to 
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Finally, to lowest order, the energy equation, Eq. (6), gives 
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Eliminating V-V between Eqs. (11) and (15), making use of Eq. (12), we obtain 
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Here, j3 = T Pq/Bq is (F times) the plasma beta calculated with the "guide-field", Bq, and 
C/3 = \//3/(l + /?)• Note that our ordering scheme does not constrain P to be either much 
less than or much greater than unity. 

Equation (15) implies that V-V ~ 0{Bq^): i.e., that the flow is a/most incompressible. 
Hence, to lowest order, we can write 



y = V0 X ^ + \4 z. 



(17) 



F. Final equations 



Let dp = Cfs dil \/\ + r, Z = hz/cp vT+T, and Vz = V^/vT+T. Neglecting the bar over 
Vz, our final set of reduced, 2-D, two-fluid, drift-MHD equations takes the form: 



^ = [^-d,Z,ij]+r^J- ^-^"^^ + ^^ V'lVz + {d,/c,) J], 
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— = [0, Z] + cp [Vz + (dp/cp) J, ^] + c| 
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+fiedp'^\U-dpY), (19) 

^ = [0>f^] - ^ {v^b,^] + [u,z] + [r,0]} + [j,^] + ^^v\u + dpTY) 

+fieV^{U-dpY), (20) 

dV 

-^ = [0, \/J + Cp [Z, iP] + /i, V V, + /Xe V'[V; + idf,/Cf3) J]. (21) 

Here, f/ = V^0, J = V^V'! ^''^^ ^ = V^Z. The four fields which are evolved in the above 
equations are the magnetic fiux-function, ip, the (normalized) perturbed z-directed magnetic 
field, Z (= bz/cjs ^/l + r), the 2;-directed guiding-center vorticity, U, and the (normalized) 



z-directed guiding-center (and ion) fluid velocity, Vz (= F-z/vT+r). The (normalized) 
2;-directed electron fluid velocity is V^ + {dp/cp) J . The quantity is the guiding-center 
stream-function. The ion stream-function takes the form 0j = (p + dpT Z, whereas the 
electron stream-function is written (pe = <P — dp Z. The above equations are "reduced" in 
the sense that they do not contain the compressible Alfven wave. However, they do contain 
the shear-Alfven wave, the magnetoacoustic wave, the whistler wave, and the kinetic-Alfven 
wave. Our equations are similar to the "four-field" equations of Hazeltine, Kotschenreuther, 
and Morrison, ^^ except that they are not limited to small values of pi. 

III. ISLAND PHYSICS 
A. Introduction 

The aim of this section is to derive expressions determining the phase- velocity and width 
of an isolated magnetic island (representing the final, nonlinear stage of a tearing instability) 
from the previously derived set of reduced, 2-D, two-fluid, drift-MHD equations. 

Consider a slab plasma which is periodic in the y-direction with periodicity length /. Let 
the system be symmetric ahont a; = 0: i.e., tp{—x,y,t) = ip{x,y,tP), Z{—x,y,t) = —Z{x,y,t), 
(p{—x,y,t) = ~(p{x,y,t), and Vz{—x,y,t) = Vz{x,y,t). Consider a quasi-static, constant--^ 
magnetic island, centered on x = 0. It is convenient to transform to the island rest-frame, 



in which d/dt ~ 0. Suppose that the island is embedded in a plasma with uniform (but 
different) y-directed ion and electron ffuid velocities. We are searching for an island solution 
in which the ion/electron fluid velocities asymptote to these uniform velocities far from the 
island separatrix. 

B. Island geometry 

In the immediate vicinity of the island, we can write 

2 

ij{x,e,t) = -'^ + ^{t)cose, (22) 

where 6 = ky, k = 27r/l, and ^{t) > is the reconnected magnetic flux (which is assumed to 
have a very weak time dependence). As is well-known, the above expression for ip describes 
a "cat's eye" magnetic island of full-width (in the x-direction) W = Aw, where w = \/^. 
The region inside the magnetic separatrix corresponds to ip > —^, the region outside the 
separatrix corresponds to ip < — !^, and the separatrix itself corresponds to ■(/' = —^. The 
island O- and X-points are located at {x,9) = (0,0), and {x,6) = (0,7r), respectively. 
It is helpful to define a flux-surface average operator: 

(/(..*.»))= /^%^? (23) 



|a;| 27r 
for ip < —^, and 

for ip > —\P. Here, s = sgn(a;), and x{s,ip,6o) = (with tt > ^o > 0). The most important 
property of this operator is that 

{[Ai'])^0, (25) 

for any field A(s, ip, 9). 



C. Island equations 

The equations governing the quasi-static island [which follow from Eqs. (18)-(21)] are: 

^ cos^ = [0-rf^Z,^] +r^5J- ^^^^^^i^ V^iV; + {dp/cp)5Ji (26) 

at Cfj 

= [0, Z] + cp [V; + {dp/cp) 6 J, ij]+clDY + /ie dp V\U -dpY), (27) 

= [0,f/] - ^ {V^[0,Z] + [U,Z] + [¥,(/>]} + [5J,iP]+fi,V\U + dpTY) 

+fi,V\U-dpY), (28) 

= [0, K] + cp [Z, ^] + fi, VV, + /ie V'[V; + {dp/cp) 6 J], (29) 

where 6J = 1 + V'^ip (the 1 represents an externally applied, inductive electric field main- 
taining the equilibrium plasma current), Y = V'^Z, U = V^0, and 

^ / 3 r \ K , , 



D. Ordering scheme 

We adopt the following ordering of terms appearing in Eqs. (26)-(29): ^ = ^^'^\ (p = 

(l)m(s,tp) + (i)'^^\s,tp,e), z = z(^\s,tij) + z(^\s,tp,e),v, = vj^\s,^ij,e), 6J = 5J^^\s,^,e). 

Moreover, V = V(°), r = r^, cp = 4°\ dp = df , /i,,e = ^lf}, «: = t^^^\ V = V^^\ 
and d\P/dt = d\P^'^ydt. Here, the superscript ^^'^ indicated an ith order quantity. This 
ordering, which is completely self-consistent, implies weak {i.e., strongly sub-Alfvenic and 
sub-magnetoacoustic) diamagnetic fiows, and very long {i.e., very much longer than the 
Alfven time) transport evolution time-scales. According to our scheme, both Z and are 
flux-surface functions, to lowest order. In other words, the lowest order electron and ion 
stream-functions, (pe = 4> — dp Z and (pi = (p + dpT Z, respectively, are fiux-surface functions. 
To lowest and next lowest orders, Eqs. (26)-(29) yield: 

^ cos^= [0(^)-rf,Z(^),^]+r^(^)5j(^)-^^2lMi±l) v^[V;(^) + (Vc,)5/^)], (31) 
at cp 



= cp [Ki^) + {dp/cp) <5J(2), ^] + cj D(2) F(i) + /i(2) dp V2(f/(^) - dp F(^)), (32) 

+fif^ V2(f/« + rf^rrW) + /x(2) V2(f/« - rf^yW), (33) 

= -M« [Vj'\,p] + c, [Z(3), ^] + ^f ) VVi^) + ^f V^[\/i^) + (dp/cp) 5J^% (34) 

where F^^) = V^Z^^), f/^^) = VV^^\ M(i)(s,^) = dcj)^^'^ / d^) , and ^(^(s,^) = dZ^^'^/di>. 
Here, we have neglected the superscripts on zeroth order quantities, for the sake of clarity. 
In the following, we shall neglect all superscripts, except for those on 0*^^^ and Z'^'^\ for ease 
of notation. 

E. Boundary conditions 

It is easily demonstrated that the y-components of the (lowest order) electron and ion 
fluid velocities (in the island rest frame) take the form Vey = x {M — dpL) and Viy = 
X {M + dpTL), respectively. Incidentally, since Vey and Viy are even functions of x, it 
follows that M{s,iIj) and L{s,iIj) are odd functions. We immediately conclude that M{s,iIj) 
and L{s,iIj) are both zero inside the island separatrix (since it is impossible to have a non- 
zero odd flux-surface function in this region). Now, we are searching for island solutions for 
which xM ^ Mq and x L ^ Lq as \x\/w -^ oo. In other words, we desire solutions which 
match to an unperturbed plasma far from the island. If V}y^ and V^ are the unperturbed 
y-directed electron and ion fluid velocities in the lab. frame, then Vjy^ — V = Mq — dp Lq and 
'^iy — V = Mq + dpT Lq, where V is the island phase- velocity in the lab. frame. It follows 
that Lq = (4°)- V;W)/rf^(l + r) and Mq = vl,%-V, where V^, = (v!^^ +rV}'J)/{l + r) 
is the unperturbed plasma E x B velocity in the lab. frame. Hence, determining the island 
phase- velocity is equivalent to determining the value of Mq. 

F. Determination of flow profiles 

Flux-surface averaging Eqs. (32) and (33), we obtain 
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(vV) + '^^^^'^ ("^^ (v^r) = 0, 

(/ij + lie) 



(35) 



and 



5' w' {W'Y) - {Y) =0, 



(36) 



where 



Cli / /X^ /^e 



. u (37) 

t(7 V -D y /Uj + Ate 

Assuming that the island is "thin" (i.e., w <^ /), we can write V^ ~ d'^ jdx^ . Hence, 
Eqs. (35) and (36) yield 

dfi (/ij r - /ie) 



M(s,V^) 



(/ij + /ie 



L(s,^) + F(s,^), 



(38) 



where 



d_ 

dip 



^ (,...(,.) ^)_(,.)i 



(39) 



and 



dip^y 'dip) 



2>(^')^)=o- 



(40) 



We can integrate Eq. (39) once to give 



*^"^^(<^*)S-<^^''^-^'°- 



(41) 



We can solve Eq. (40), subject to the constraints that F be continuous, F = inside the 
separatrix, and F —^ s Fq as |x|/w — > oo, to give 

rV- dip I /■- °° dip 



Fis,iP) = sFo 



-^ (a;4) / 7-^ (a;4) 

outside the separatrix. Note that x F ^ \x\ Fq as \x\/w — > oo. 

In order to solve Eq. (41), we write ip = —ip/^, ((■ ■ ■)) = (■ ■ ■)w, X = x/w, and L 
L/{Lo/w). It follows that 



(42) 
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Suppose that 6 <^1. In this case, L{s,ip) takes the value s/((X^)) in the region outside the 
magnetic separatrix, apart from a thin boundary layer on the separatrix itself of width 6 w. 
In this layer, the function L{s, ip) makes a smooth transition from its exterior value (which 
is stt/4 immediately outside the separatrix) to its interior value 0. We can write 

L{sJ) = s(^j^^+l{y)y (44) 

where y = {ip — 1)/S. It follows that 

dH 3 ^ 

since {{X^))^-^ = 4/7r, and ((X'^))j;=i = 32/3 tt. Hence, the continuous solution to Eq. (41) 
which satisfies the appropriate boundary conditions is 



L{sJ) 



1 n / 3 ^-1 

— exp 



(46) 



((X2)) 4 -^ V8 5 

in the region outside the separatrix {i.e., ip > 1). Of course, L{s,ijj) = in the region inside 
the separatrix {i.e., ip < 1). 

G. Determination of island phase-velocity 

Let 6 J = 6Jc + SJs, where 6Jc has the symmetry of cos^, whereas 6Js has the symmetry 
of sin^. Now, it is easily demonstrated that 

{5.J, sine) = ^{x[5.J,^]). (47) 

Hence, it follows from Eq. (33) and (38) that 

(«-.)^-(-±/^A((.=)|£-.(.nf-W.). (48) 

Now, for an isolated magnetic island which is not interacting electromagnetically with any 
external structure, such as a resistive wall, the net electromagnetic force acting on the island 
must be zero. This constraint translates to the well-known requirement that 

12 
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{6JsSme)d^ = 0. (49) 



Using Eq. (48), this requirement reduces to the condition 

hm I (x^) —— - 2 (x^) -— - (x)F I oc hm 

x/w^oo \ dlp^ dip J x/w^oo 



2 d fl d{xFy 
dx\x dx 



= -Fo = 0, (50) 

since xF -^ \x\Fq as \x\/w -^ oo. Hence, we conclude that Fq = [i.e., F{iIj) = 0, 
everywhere] for an isolated magnetic island. 
It follows from Eq. (38) that 

M(.,V^) = - ^^/^-^~f^ L(g,^). (51) 

Hence, Mq = -[dpi^jiiT - fie)/ifJ'i + /^e)] -^o- Recalling that Mq = V^^y - V, dpL^ = 
iV,f - V;T)/(1 + ^)' V^ = ^Ely + dpTL,, and I/W = \4°i, - d^L,, we obtain the 
following expression for the island phase-velocity: 

V = ^»^y +l^eVey ^^2) 

In other words, the island phase- velocity is the viscosity weighted mean of the unperturbed 
ion and electron fluid velocities. Hence, if the ions are far more viscous then the electrons, 
then the island propagates with the ion fluid. In this case, the ion fluid velocity profile 
remains largely unaffected by the island, but the electron fluid velocity profile is highly 
sheared just outside the island separatrix. The opposite is true if the electrons are far more 
viscous than the ions. This is illustrated in Fig. 1. 
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FIG. 1. Velocity profiles as functions of x, at constant 9, evaluated on a line passing through 
the island 0-point (i.e., at 9 = 0) in the island rest frame. The 0-point lies at x = 0. The island 
separatrix is indicated by a vertical dotted line. The solid curves show the normalized ion fluid 
velocity profile: [Viy — V)/{V^ — V}y ). The short-dashed curves show the normalized electron 
fluid velocity profile: {Vey — ^)/(^„ ~ ^ey ). The long-dashed curves show the normalized E x B 
velocity profile: (VEBy — ^)/(^V ~ ^ey )■ ^^^ left-hand panel shows the case of viscous ions: 
fie/lJ-i = 0.1, r = 1., and 6 = 0.2. The right-hand panel shows the case of viscous electrons: 
fJ'i/fJ'e = 0.1, r = 1., and 6 = 0.2. 

We have now fully specified the ion and electron stream- functions, 0j and (f)e, respectively, 
in the island rest frame. In fact, 0j = inside the separatrix, and 



d(l>i{s,i[j) ^ (0) _ (0). fij L{s,ij) 

J„l. K*^ iy ^ ey J , ,, 



(53) 



outside the separatrix, where the function L{s,iIj) is specified in Eq. (46). Likewise, 0e is 
zero inside the separatrix, and 



dip *^ "^^ Hi + fXe W 



(54) 



outside the separatrix. Note that the stream-functions and their first derivatives are every- 
where continuous, which implies that the ion and electron fluid velocities are everywhere 
continuous. 
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H. Determination of ion polarization correction 

It follows from Eq. (33) that 

where I{s,iIj) is as yet undetermined. The function Hli)) is zero for i) < 0, and unity for 
■(9 > 0. Here, we have made use of the fact that outside the separatrix L{s,tp) ~ sLq/{x^), 
and M{s, ■?/') ~ s Mq/{x'^), apart from a thin boundary layer on the separatrix itself. It turns 
out that we do not need to resolve this boundary layer in order to calculate the total ion 
polarization current. However, we do have to include the net current flowing in this layer 
in our calculation of the total current. ^^'^^ Flux-surface averaging Eqs. (31) and (34), we 
obtain 

e'w' (V^^Je) - ((5Je) = -v~' ^ (cos^), (56) 

where 



"j / f^i f^e (K7\ 



W y/r] ]l fli + He 



Equation (56) implies that 



1 d^ 

(5J,)^^-i — (cos^), (58) 

apart from in a thin boundary layer on the separatrix of width ew. Here, we are assuming 
that e <^ 1. It is easily demonstrated that the deviation of {6Jc) in the boundary layer 
from the value given in Eq. (58) makes a negligible contribution to the total ion polarization 
current. Hence, we shall treat Eq. (58) as if it applied everywhere. 
Equations (55) and (58) give 

Note that this current profile contains no discontinuities or singularities. 
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The island width evolution equation is obtained by asymptotic matching to the region 
far from the island.^ In fact, 



-oo 



A'^ = -4 {5. J, cosO) dij, (60) 

where A' is the conventional tearing stability index. ^^ It follows from Eqs. (59) and (60) 
that 

Ailv^^"'^^ ((!)) j#l((XTJ "^ 

8dwr-{(s^^_ (61) 

1] dt J-i ((!)) ^ ^ ^ 

Performing the flux-surface integrals, whose values are well-known, ^^ we obtain the following 

island width evolution equation: 

it) iv -^°^ 

ri dt ' ^'^^ {W/Af 



0.823 dW ,, {V - \4'^ J {V - V,\ , 

A' + 1.38^ .:.,... ^^- (62) 



Here, W = A:W vs, the full island width. The ion polarization current term (the second term 
on the r.h.s.) is stabilizing when the island phase- velocity, V ^ lies between the unperturbed 
local E ^V velocity, V^^y, and the unperturbed local velocity of the ion fluid, V^y }^ 

IV. SUMMARY AND DISSCUSION 

A set of reduced, 2-D, two-fluid, drift-MHD equations is developed. This set of equations 
takes into account both electron and ion diamagnetism (including the contribution of the 
ion gyroviscous tensor), as well as the Hall effect and parallel electron compressibility, but 
neglects electron inertia and the compressible Alfven wave. For the sake of simplicity, the 
plasma density is assumed to be uniform, and the ion and electron temperatures constant 
multiples of one another. However, these constraints could easily be relaxed. 

Using our equations, we have derived a complete and self-consistent solution for an iso- 
lated magnetic island propagating through a slab plasma with uniform but different ion and 
electron fluid velocities. Our solution is valid provided that the ordering scheme described 
in Sect. HID holds good, and the island width W is sufficiently large that 
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^^d^J^^ (63) 



[i.e., 6 <^ 1), and 



^^d,[J^ (64) 

(i.e., e< 1). 

Note that the ordering scheme described in Sect. HID imphes that cj* ^ k\\Cs^k\\VA 
where cu* is a typical diamagnetic frequency, Cg the sound speed, and va the shear-Alfven 
speed. Here, k\\ must be evaluated at the edge of the island. This scheme differs from that 
adopted in Ref. 16, for which k\\ Cs ^ uj^f. It turns out that our ordering scheme permits a 
much less complicated calculation of the flow profiles around the island than that described 
in Ref. 16. 

Within our solution, the ion and electron fiuid velocity profiles are uniquely determined 
in the vicinity of the island [see Fig. 1]. These profiles are everywhere continuous and 
asymptote to the unperturbed fiuid velocities far from the island. Incidentally, the inclusion 
of electron viscosity in both the Ohm's law and the plasma equation of motion is key to the 
determination of continuous velocity profiles. ^^ 

The island phase-velocity is uniquely specified by the condition that there be zero net 
electromagnetic force acting on the island [see Eq. (52)]. It turns out that the phase- velocity 
is the viscosity weighted mean of the unperturbed ion and electron fiuid velocities. In this pa- 
per, we have adopted phenomenological diffusive ion and electron viscosity operators, which 
are supposed to represent anomalous perpendicular momentum transport due to small-scale 
plasma turbulence. 

The ion polarization current correction to the Rutherford island width evolution equa- 
tion is found to be stabilizing when the island phase-velocity lies between the unperturbed 
ion fiuid velocity and the unperturbed E x B velocity [see Eq. (62)].^^ It follows, from our 
result for the island phase- velocity, that the polarization term is stabilizing when the anoma- 
lous perpendicular ion viscosity significantly exceeds the anomalous perpendicular electron 
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viscosity [see Fig. 1, left panel]. Conversely, the polarization term is destabilizing when 
the electron viscosity significantly exceeds the ion viscosity [see Fig. 1, right panel]. ^^ Note, 
however, that in order for the electron viscosity to exceed the ion viscosity, the electron 
momentum confinement time would need to be at least a mass ratio smaller than the ion 
momentum confinement time, which does not seem very probable. Hence, we conclude that 
under normal circumstances the polarization term is stabilizing. 
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